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case allowing an exact time-optimal realization of environment-mediated non-local system unitaries. 
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I. INTRODUCTION 

The discussion of optimal control protocols for evolving 
quantum systems and the relevance of transition speeds 
and their limits can be traced back to the dawn of quan¬ 
tum theory as can be seen in the example of Fermi’s 
golden rule for the mean transition rate into an orthogo¬ 
nal energy eigenstate. While touchiM fundamental prob¬ 
lems about computation per se [iHsI, it is also of tech¬ 
nological interest, since minimization of computing time 
amounts to minimizing the interaction time with detri¬ 
mental environmental noise sources. The speed optimiza¬ 
tion can refer to the transformation of a given initial to 
a target state, in particular the evolution into an (arbi¬ 
trary) orthogonal state, or it can refer to the implementa¬ 
tion of a desired target unitary, which is supposed to act 
on an arbitrary (unknown) quantum state. The existence 
of a physical solution to the optimization problem and the 
formulation of finite speed limits generally require a set 
of preset constraints which represent practical limitations 
encountered. They can be of energetic nature restricting 
the size of the Hamiltonians in some way or refer to the 
type (subspace) of Hamiltonians that one can implement. 
The latter kind of constraints can be added by Lagrange 
multipliers giving rise to quantum brachistochrone 
curves as optimal trajectories of states or unitaries. 

Besides of this, a fundamental constraint addresses ex¬ 
ternal background fields that cannot be manipulated and 
give rise to a free (natural) Hamiltonian Hq, which must 
be taken as given. Even in the absence of restrictions 
on the implementable control Hamiltonians Hq {t) (which 
would suggest to simply substitute Hq = He — Hq), the 
mentioned energetic constraints limit the size of imple¬ 
mentable control fields, hence the generic Hq cannot be 
ignored. This is the motivation for the quantum Zermelo 
navigation problem which, inspired by a result in Q , was 
formulated in Q and recently solved in [1, @ . 

Roughly speaking, external forces as caused by a back¬ 
ground field are incorporated into the problem descrip¬ 
tion as a geometric property by defining an appropriate 


measure of distance with respect to which one then freely 
takes a shortest path. In the example of the Zermelo nav¬ 
igation problem, the solution consists in following the 
geodesics of a Randers metric, which is a special type of 
Finsler metric and can be derived from Riemann metrics 
by adding a linear term. We refer to [1, @] and references 
therein for a description of the mathematical background. 

We here consider the navigation of unitary operators. 
Starting from a given initial unitary Ui, the task is to find 
a time-dependent Hamiltonian H(t) that implements a 
desired final unitary t/p in shortest possible time. The 
total Hamiltonian 

H{t)=Ho+Hcit) (1) 

is assumed to consist of a constant part Hq that cannot 
be altered and a time-dependent part Hc{t) which can 
be controlled by the experimenter without limits except 
that the size of the control is bounded, so that we set a 
fixed Tr[ILQ(t)]. 

The reference to Zermelo originates in the classical 
problem of finding a time-dependent heading of a ship 
or aircraft that, starting from some initial locatiorq navi¬ 
gates it to a destination in shortest possible time @. Hq 
here corresponds to a wind or current, while a constant 
Ti[HQ{t)] corresponds to the assumption of “full speed 
ahead” at any time. In the classical problem, the opti¬ 
mal route under constant wind or current is a straight 
line from start to finish, which is achieved by keeping a 
constant correction angle of the vessel’s heading with re¬ 
spect to its destination to compensate the drift off. More 
general versions of the classical Zermelo navigation prob¬ 
lem that include obstacles and a wind depending on loca¬ 
tion and time have also been of recent interest in animal 
behavior research 0 and computer science 0. 

Contrary to intuition, the solution of the correspond¬ 
ing quantum problem is not a constant Hc(t) = Hc which 
would give rise to a (single-exponential) one-parameter 
subgroup of unitaries U{t) = Te~^ —Q-i{Ho+Hc)t^ 

but instead an explicit time-dependence such that H be- 
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comes time-independent in the interaction picture 

-ffcW =e-‘^“‘i?c(0)e‘^«‘, (2) 

and for which the time-optimal curve of unitaries 
disentangles according to U(t) = = 

g-iifotg-iffc(o)t _ exponentials. 

(Here, T denotes positive time ordering, and we have set 
h=l.) Time-independence as in the classical case is thus 
only recovered for [Hq, .Hc(O)] = 0. 

Zermelo navigation of quantum states has been solved 
in [l^ also deriving a limit on the passage time, by trac¬ 
ing it back to the navigation of unitaries. The energetic 
constraint here takes the form of a fixed variance (AHq). 
The special case of time-independent Hq has been dis¬ 
cussed in [I^, placing focus on two-level systems and 
demonstrating that the treatment is astonishingly non¬ 
trivial even with these restrictions. Interestingly, Eq.(5) 
in [I^ transforms the state navigation problem I'I'f) 
= [/(t)l'I'i) (in its general form) into that of hitting a 
(counter)moving target, = e'^‘’‘|'I'F), which 

suggests a close connection to pursuit problems: it ad¬ 
vises to move “straightly” to where the target will be 
at the time of arrival. This is analogous to the classical 
case of a dog who should straightly move towards the ex¬ 
trapolated target (her walking master’s) location, rather 
than intuitively following her snout’s direction heading 
towards the instantaneous location of the target result¬ 
ing in a longer so-called “Hundekurve”. 

Referring to time-independent systems with Hamilto¬ 
nian H, |14l | relates speed limits for the implementation 
of unitaries to a family of bounds on orthogonality times, 
which are all of the canonical form t± > ^, includ¬ 
ing the Margolus-Levitin theorem [I^ {E = {tp\H\ip) — 
Eo, Eo being the lowest eigenvalue of H), while relation 
to the time-energy uncertainty relation (Mandelstam- 

Tamm bound [l^ , E = ^— (ip\H\iljy) is also dis¬ 
cussed. 

Here, we show how the Zermelo navigation ([2]) can be 
applied to a case of non-Markovian evolution of open 
quantum systems such that it allows an exact treat¬ 
ment of a set of quantum systems coupled to a com¬ 
mon environment. In contrast to the original motivation 
for realizing a time-optimal evolution, namely to mini¬ 
mize the temporal accumulation of environment-induced 
decoherence effects, we demonstrate that Zermelo time- 
optimality can enable a deliberate environmental cou¬ 
pling, which here serves to mediate a desired system 
coupling. This intended utilization of non-Markovian 
(memory) properties of the environment thus contrasts 
the case, where environmental interactions are accom¬ 
panied with irreversible decoherence. Bath-induced sys¬ 
tem interactions have been discussed for several years, 
cf. e.g. [3, [l^. The main focus of these Markovian ap¬ 
proaches has been on the preparation of entangled many- 
body steady states [20l - [^ for dissipative quantum com¬ 
putation and continuous quantum repeaters [2^. 


In contrast, here we consider a navigation of unitaries 
acting on a set of systems and a shared bath. Start¬ 
ing from the identity, Ui = /, the navigation periodically 
passes unitaries U{tm) in which system and bath parts 
are refactorized. Although the systems never interact 
directly, the system part of the target unitary contains 
an environment-induced pairwise system-coupling which 
is decoherence-free. All results are independent of the 
system and bath states and apply exact analytical ex¬ 
pressions which do not rely on assumptions such as the 
weak coupling perturbative approach, rotating wave or 
Markovian approximations. 

The present work builds on [26l - l^ . who demonstrate 
the preparation of nonclassical states such as the GHZ- 
state, cf. also [l^. However, the scheme presented there 
requires setting Hq = 0 to allow for an analytic solution. 
Below we will show that (a) it is remarkably the Zermelo 
navigation m that enables us to overcome this limita¬ 
tion, allowing for arbitrary Hq. As a consequence, the 
described protocol of generating bath-mediated system¬ 
coupling is tim e-op timal, a result thus also holding for 
the schemes in |26l - l^ as a special case of vanishing Hq. 
Based on this, we construct (b) a state-independent speed 
limit for the refactorization time. Finally, we address 
the question of reaching an arbitrarily given target U-p- 
While a single Zermelo navigation U (t) allows in our case 
merely the creation of bath-induced pairwise unitary sys¬ 
tem interactions, the fact that U{t) is a product of two 
non-commuting exponentials suggests the reachability of 
a given Up by concatenation as realized by Hamiltonian 
resets. The two-exponential form of U(t) thus halves the 
number of required resets. We will (c) provide a numer¬ 
ical proof of principle demonstration, where the systems 
are three two-level systems (“qubits”), and Up is a Toffoli 
gate. Apart from being optimal for quantum error correc¬ 
tion, the Toffoli gate alone is universal for reversible com¬ 
puting and, together with the Hadamard gate, forms a 
universal set of gates for quantum computing [30l| . While 
the three qubit-setup here serves as simplest case going 
beyond pairwise system interactions, it is hence already 
sufficient as building block for universal quantum com¬ 
putation. 

This work is organized as follows. After this introduc¬ 
tion, the realization of an exact decoherence-free bath- 
mediated unitary system coupling by means of Zermelo 
navigation is presented in Sec. [iTl Sec. Ill Al defines the 
general setting and conditions. Sec. IHBI introduces the 
example of a linear coupling to a bosonic bath, for which 
the speed limit is derived in Sec. Ill Cl A numerical ap¬ 
proach to the problem of navigating to a general target 
is given in Sec. lHIl As a setup relevant for quantum com¬ 
puting, the special case of three qubits is considered for a 
Toffoli gate in Sec. IHI Al and for a repeater relay station 
in Sec. IHIBI An approach to the navigation to a general 
target that does not rely on the requirements of the pre¬ 
vious sections is outlined in Sec. IIVI Finally, a summary 
and outlook on future work is provided in Sec. El 
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II. BATH-INDUCED UNITARIES 

A. Zermelo navigation for bath-coupled systems 

Consider N systems j = 1,..., N and one bath B as 
depicted in Fig. [T] with factorized individual system-bath 


ordered product according to (sTj 

= ^1 (^2, ti )Te~‘'^**i^ 

= Ui{t2 


(7) 

( 8 ) 

(9) 
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FIG. 1: A set of systems Sj interacting with a shared bath 
(environment) B. An extension with modulated system-bath 
interaction strengths, where each term in is multiplied 
with a will be discussed in Sec. IIIII below. 


interactions, 

H = ^s+^B+^SB, (3) 

N N 

Hs = Y. -^SB = E ® ’ (4) 

i=i i=i 

where Sj are Hermitian operators of system j (hence 
[Sj, 5fc] = 0) and Bj are Hermitian operators of the bath 
(hence [Bj,Hfe] 0). 

The division in the quantum Zermelo navigation prob¬ 
lem of the Hamiltonian (HD into an uncontrollable “wind” 
part Hq and a controllable part Hq is not fixed and refers 
to the operators one can and wants to alter. In our case 
m of a composite quantum system, a natural assump¬ 
tion is that one can modify operators of the system but 
not those of the bath. We may thus identify in ([3|) 

Ho = Hs, Hc = Hb + Hsb, (5) 

for which ([2]) with Hc{^) =Hc becomes 

N 

( 6 ) 

i=i 

where Sj{t) = Uj{t)SjUj{t) with Uj{t) =e“'^^‘. We see 
that only the individual system operators Sj have to be 
navigated in their respective local winds Hj , whereas all 
bath operators Hb and Bj remain unaltered. 

Given two general non-commuting and explicitly 
time-dependent Hamiltonians Hi{t) and H 2 {t), defin¬ 
ing lJi{t 2 ,ti) = Te“‘7ti using the conven¬ 

tion Cl(< 1 ,^ 2 ) = Ul{t 2 ,ti), we can disentangle the time- 


where the auxiliary time can be chosen as desired. 
Applying ([5]) with ti =0, t 2 = t, Hi = Hs + Hb, 772(t) 
= 5i(T) G) Bj gives for the total evolution operator 

U{t) = U{t,0) 

U{t) = Te"'7od^ff5) 

= UsmBmsBit), (10) 

UsBit) = Te-'7odrAsB(r) (H) 

where Us(t) = e“'^®‘, C/b( 7) = and HsBiT) = 

5j <8 Cj(r) with Bj^r) = ul{T)BjUBiT). [In phys¬ 
ical terms, this dynamic time dependence of 77 sb(t) via 
Bjij) refers to the Dirac interaction picture and must 
be distinguished from the explicit time dependence of 
Hc{t) via Sj{t) introduced in dH]). The latter refers to 
the Schrodinger picture and describes the actual naviga¬ 
tion carried out by the experimenter.] 

Eq. (ITT]) defines some Hermitian operator HeS, which 
can be expressed by means of the Magnus series expan¬ 
sion [s^. iJeff itself is time-dependent, and we have left 
a factor t for dimensional reasons. Since 

N 

[HsB{ti),HsB{t2)] = Y, Sj[Bj{ti),Bk{t2)]Sk, ( 12 ) 

j,k=l 

in the special case when [Bj(ti), I3k{t2)] are c-number 
functions, only the first two terms of the Magnus series 
are nonzero, and we obtain a closed expression 

1 /■* 

Hes{t) = - dTHsBir) 

Jo 

— ^ J dti J dt2[HsB{ti), HsBih)] (13) 

N N 

= Y SjK:jk{t)Sk, (14) 

i=i 

where 

■C (l) = J dTCj(T), (15) 

Kjkit) = [ dti [ dt2[Bj{ti),Bkit2)]. (16) 

Jo Jo 


B. Linear coupling to a bosonic bath 

A simple example in which the Magnus series breaks 
off at second order is a bath of bosonic modes [26l - [^ . 
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The bath parts of Hsb in Q 

Bj = ^ ^ -Ajibi + h.a. (17) 

i 

are linear combinations of the bosonic mode operators 
bi and h.a. denotes the Hermitian conjugate. A bath 
Hamiltonian 

Hb = ^ (18) 

i 

then gives the time evolution Bj{t) = + 

h.a. Using [S/,SJ,] = Sw we can verify that the above- 
mentioned c-number condition is fulfilled, [I3j{ti), Bkit 2 )] 
= 2ilm^; If the Aji couple only 

to harmonics of some given frequency w, i.e., they are 
nonzero only for uji = Ilo, I € N, then for discrete times 

tm = mAt, At =— (19) 

(jJ 

(to e N) the integral (ITS!) vanishes, Fj{tm) = 0, whereas 
(HSI reduces to Kjk = Kjkitm) = 

Combining the Aji to a matrix A and defining a diagonal 
matrix Jl with elements U/;/ = —6iii/{loj) we can write Kjk 
= Re(Ar2Al)jfc and combine the Kjk to a matrix k. 

In summary, the total evolution operator (HOI) factor¬ 
izes at the discrete times tm according to 

tj{mAt) = Q-im^tHs^-iraAtS n-S ^ ^ ^20) 

where S ■ k ■ S = 'Y^^k=i SjKjkSk- The first two expo¬ 
nentials are a result of the Zermelo navigation, whereas 
the last one describes the free evolution of the decoupled 
bath. 


C. Speed limit 

In our scenario, a set of N systems interact with a com¬ 
mon bath such that the overall unitary evolution factor¬ 
izes for discrete times m into system and bath parts. 
The bath nevertheless affects the system evolution by 
means of the exponential containing the^ in ( 1 ^ 01 ) . This 
is in contrast to Markovian bath effects [l3] , where infor¬ 
mation dissipates irreversibly into the environment thus 
preventing the realization of a unitary transformation. 
Unlike the system evolution, the evolution of the bath it¬ 
self in ( 1 ^ 01 ) is not affected. The Zermelo navigation serves 
to remove any time dependence from the control Hamilto¬ 
nian in the interaction picture which allows a straightfor¬ 
ward evaluation of the obtained target unitary by means 
of the Magnus identity. Due to the need of a well-defined 
coupling, the bath here acts as an ancillary system such 
as a field confined by a cavity rather than an uncontrolled 
environment. 

Obviously, Tr[HQ{t)] =Tri... TrArTrB[i7c(^)] i® unde¬ 
termined but fixed because it is time-independent. In 


contrast to the finite-dimensional systems for which the 
quantum Zermelo navigation problem has been discussed 
so far, the bath considered here is an infinite set of 
infinite-dimensional systems. The Zermelo solution of 
the problem does not contain any reference to the sys¬ 
tem size, however, which suggests that it holds in our 
more general setting too, cf. also a comment in [^. (We 
may truncate the bath with respect to the number of 
relevant bath modes and at some sufficiently high occu¬ 
pation number for each of these modes without changing 
any physical properties of the model.) In (1^ . the in¬ 
finite bath is decoupled from the unitary evolution of 
the N systems, where the latter may be finite or infi¬ 
nite dimensional. Restricting the modulation (imposed 
time-dependency) to the systems renders the method in¬ 
dependent on any need to manipulate the bath. 

Since we apply Zermelo navigation, the described 
method of bath-induced realization of the above type 
of unitaries is time-optimal with respect to the local 
“winds” Hj. Substituting a; = ^, where Eq is the ground 
state energy of a harmonic oscillator [i.e., if we set Hb = 
fej(h-l-4)], the minimum time required for refactorization 
of the unitary into system and bath parts can be written 
in the canonical form mentioned in the beginning 

nh 

tmin = (21) 

Eo 

A single harmonic oscillator with fundamental frequency 
oj here serves as a minimum system suitable to act as 
bath in the described way. This minimal implementa¬ 
tion time is finite but small for large Eq. An observer 
on a large time scale, on which she cannot resolve tniim 
may thus be unaware of it and even the existence of a 
background bath. Such an observer would simply wit¬ 
ness a time-continuous evolution of the systems, whose 
bath-mediated interactions would exhibit as direct pair¬ 
wise system interactions. 

III. CONCATENATIONS 

So far we have i gn ored the problem of reaching a given 
target unitary. In Q, a method is presented to calculate 
the initial He (0) such that U (t) after some time t reaches 
the target Up, U{t) = UpUi, so that U{0)Ui = IJi and 
U{t)Ui = Up. In our context, we face the additional re¬ 
striction to pairwise couplings S-k-S in (1201) . The system 
part of (1201) is itself a product of two factors, Usystemirnr) 
= by i = SjKjfSf and B 

= Hj, respectively, where r is fixed by (fTOl) . This 
restricts the set of implementable A^-system target uni¬ 
taries Usystemitm = bt,t) = Up, even if the Sj can be cho¬ 
sen at will. To generalize the set of reachable Up, we may 
concatenate individual Zermelo navigations by adjusting 
at times tm = ''’fc [by means of u), cf. (HH)] t to 

Tm+i, Kjj' to [by multiplying each term in iJsB 
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in o with a system-bath interaction strength 
> 0], and reset the Sj{t) to Sj. This gives Usystemitm) 
= with Ai!^^ = ||ifcj|-iifc, Ofc = 

llAlkfe, Pk = Tk, and given Ak = with 

li'y\Tk / T)Kjj'. Here, we factored out ||4p = 
Tr(H’l’H) for convenience, and without loss of generality 
we assume that ||B|| = 1. The special case and 

hence ~ Kjji generates an alternate product for which 
A^^'^ =A^-^\ which can approximate any UF = e~'^, with 
H being a member of the algebra spanned by the multi¬ 
commutators from A and B (33| , where the ak and Pk can 
be determined from Up [s^. The resulting concatenated 
navigation represents a (piecewise differentiable) evolu¬ 
tion from I to Up under the constraint of fixed A^^'> 
(or A^^'‘) and B. The Zermelo-optimality holds only for 
the segments ', during which the two fac¬ 
tors and are hence implemented at once 

within Tfe rather than in consecution, sparing the need to 
switch between the generators Ak and B. 

The number N of systems can be arbitrary including 
infinite. In the fundamental case iV = 1, in which no 
division into subsystems is assumed, the effect of the bath 
reduces to “dressing” the system evolution, Ak^S'^. For 
iV > 1, the individual systems j and their bath-coupling 
as defined by Hj and Sj can be set arbitrarily and do 
not need to be the same for each j. For N > 2, one 
may argue that even if we can choose the Sj at will, 
they are local operators as are the Hj. However, since 
A (or Ak) is 2 -local (i.e., it represents a pairwise system- 
coupling), and each commutator of a fc-local operator 
with A generally yields a (fc-|- l)-local operator unless the 
number of systems has been reached, the concatenation 
builds up A^-system unitaries, analogous to a quantum 
circuit consisting of one- and two-qubit gates. 

A. Quantum circuits with qubits 

A “continuous variable” example is a set of har¬ 
monic oscillators with parametric-amplifier type cou¬ 
plings, Sj = hj, for which the bath induces Kerr and 
cross-Kerr nonlinearities SjSk = hj ® hk- Of special in¬ 
terest for quantum information processing are two-level 
systems, however, for which we now provide some explicit 
examples. If ihg and He are given appropriately. Up may 
represent a desired gate without need of concatenations. 
An example for iV = 2, where Hq and He commute, is the 

CNOT-gate [H Hcnot = 6 "*^e‘f 

which can be implemented if Si = —Hi = d'^\ S 2 = —H 2 
= ^^i \ and ki 2 = \. 

If the Hj and Sj are fixed by the experimental setup, 
concatenations are required to implement a desired Up. 
As a concrete example, we consider three qubits with ran¬ 


domly given Hj and Sj , whose interaction strengths with 
a shared bosonic bath can be altered [by multiplying each 
term in ihsB in (|4|) with a factor Uj ^ 0 over time inter¬ 
val Tk as mentioned above], and where the bosonic bath 
is eliminated by Zermelo-navigation of the Sj leading to 

products Hsystem(<m) =nr=i for the 

A-system unitary as described. The protocols are thus 
sequences of these finite time periods r^, over which the 
system-bath interaction strengths of the qubits j are mul¬ 
tiplied by the respective . Specifically, we consider (a) 
a synchronous protocol, generating an alter¬ 

nate product as mentioned above. An alternative is (b) 
an asynchronous protocol consisting of concatenations of 
8 -step cycles, where all subsets of the qubits are brought 
in contact with the bath as follows: over ri, all qubits 
are detached from the bath, =0, after which qubit 1 , 
then qubit 2, and then qubit 3 alone is brought in con¬ 
tact with the bath over r2, T3, r4, respectively. After this, 
only qubits 1 and 2, then only qubits 1 and 3, and then 
only qubits 2 and 3 are brought simultaneously in contact 
with the bath over rs, tq, T 7 , respectively. Finally, all 3 
qubits are brought in contact with the bath over rg. After 
this, at time tg = X]fc=i the cycle is repeated from the 
beginning over different times rg • • ■ rig, and with differ¬ 
ent interaction strengths for those qubits in bath contact, 
and so forth, until a maximum number n of exponential 
factors in C/gystem has been reached. The qubit-bath in¬ 
teraction strengths of those qubits j which are in bath 
contact over Tk are adjusted synchronously, iJ,p = > 

0 as in (a). In both protocols (a) and (b), we minimize 
the (squared) operator distance D= ||17system — 17 f||^ 
of the implemented unitary lAystem and a given target 
Up by gradient descent. The minimum distance reached 
for a given number of exponential factors n is shown in 
Fig. [2] for the case where Up is either the CNOT (see 
Sec. IIIIBl below) or the Toffoli (i.e., C^NOT) gate [13 
for both types of protocols. The choice of a Toffoli gate 
is motivated by its relevance for quantum computing, 
but analogous results can be obtained for a random gate. 
Fig. 121 demonstrates that a desired target Up can gen¬ 
erally be reached for both types of protocols as soon as 
n surpasses a threshold [> (2^)^ = 64, albeit smaller n 
may suffice for specifically given flj, Sj, and Up as in 
the closed expression for Hcnot mentioned above]. An 
example for each type of protocol is illustrated in Fig. |31 
The choice of a synchronous protocol for the Toffoli gate 
and an asynchronous protocol for the CNOT gate is here 
irrelevant, since each gate can be generated with both 
types of protocols. 

The asynchronous protocol requires a system-resolved 
switching of the coupling but exhibits a more robust con¬ 
vergence. This is illustrated in Fig. |4| where in contrast 
to Fig. 121 the three qubits and their couplings are identi¬ 
cal, i.e., there is one single (randomly given) H, so that 
the local winds (acting in the three-qubit Hilbert space) 
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FIG. 2: Minimized squared operator distance D = ||[/system — 
?7p||^ as a function of the number n of exponential factors 
in [/system for concatenated Zermelo navigations with three 
qubits j = 1,2,3 sharing a bosonic bath. The plots show 
the implementation of a [/p = CNOT gate between qubits 1 
and 3 as well as a [/p = Toffoli gate, using the synchronous 
and the asynchronous protocol. The CNOT gate refers to 
the quantum repeater relay station depicted in Fig. The 
local winds Hj and couplings Sj are randomly given for each 
individual j. 



FIG. 3: Illustration of example protocols consisting of con¬ 
catenated Zermelo navigations of qubits 1,2,3. Shown are the 
strengths ||j4j;|| of the system-bath interaction (the navigator’s 
“engine power”) over the time intervals rt for those qubits in 
bath contact during Tk- Upper plot: synchronous protocol 
generating a Toffoli gate (n = 92), lower plot: asynchronous 
protocol generating a CNOT gate between qubit 1 and 3 {n 
= 84). The CNOT gate refers to the quantum repeater relay 
station shown in Fig. 


are Hi = H ® I ® I, H2 = I ^ H 0 I, H3 = I ^ I H, 
and the local couplings Sj are analogously given by one 
single (randomly given) S. Fig. 0] demonstrates that the 
synchronous protocol fails in this case, in contrast to the 
asynchronous protocol. 


FIG. 4: Same analysis as in Fig. [2] except that the randomly 
chosen local winds Hj and couplings Sj are identical for each 

j- 


B. Quantum repeater relay station 


Depending on the application, the system coupling 
may follow a restricted topology. For example, the sys¬ 
tems may be arranged in a chain, and bath-induced cou¬ 
pling is possible only for nearest neighbors. In the sim¬ 
plest case, three qubits form a linear array as shown in 
Fig. in with a presence of direct couplings K 12 and K 23 , 
but absence of an end-to-end coupling, K 13 = 0. This 



^12 



^23 


CNOT 13 (ki 3 = 0 ) 



FIG. 5: Quantum repeater relay station: qubits 1 and 3 are 
coupled to qubit 2 but not to each other directly. The inter¬ 
mediate qubit 2 can be eliminated by concatenated Zermelo 
navigations to generate a direct coupling between qubits 1 and 
3. An example protocol for generating a CNOT gate between 
these qubits is shown in the lower plot of Fig. [3] Analogous 
results can be obtained for a random gate. 


setup resembles our original problem with the role of the 
bath adopted by the intermediate qubit 2. Similar to 
the elimination of the bath by Zermelo navigation, the 
intermediate qubit 2 can be eliminated by concatenated 
Zermelo navigations to generate a (qubit 2 - mediated) 
direct coupling between qubits 1 and 3 alone and thus a 
desired target gate Up . Since nesting of this procedure 
on larger time scales for a successive elimination of blocks 
of subsystems is (aside from the open system dynamics) 
reminiscent of quantum repeaters (SSH^ , we refer to the 
setup depicted in Fig. [5] as a quantum repeater relay sta¬ 
tion. 
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IV. GENERALIZATIONS 

A. Motivation 

Although concatenating Zermelo navigations as dis¬ 
cussed in Sec. Illll is a straightforward option for reaching 
a target, similar to the construction of a quantum circuit 
from individual gates, the concatenations as a whole are 
no longer time-optimal, since they introduce refactoriza¬ 
tions of the overall system-bath unitary at the interme¬ 
diate times tm = Ym^=i '’’fe artifact of this approach. 
The optimal system-bath trajectory would be a brachis¬ 
tochrone U (t) that starts at the identity and refactorizes 
only at the time of arrival. Even if we choose (concate¬ 
nated) Zermelo navigations, we must ensure (a) the ex¬ 
perimental realization of the required modulations Sj (t) 
(and their resets), (b) the presence of a linear coupling 
to a bosonic bath with the desired exclusive coupling to 
harmonics of a fundamental frequency w (and its adjust¬ 
ment), and (c) detailed knowledge of all operators Hj, 
Sj , and coupling coefficients Kjji , along with an absence 
of operator noise and decoherence. The realization of 
these requirements can be challenging. 

To avoid these problems, we now consider an alter¬ 
native that relies on measurements alone, i.e., a “closed 
loop” - scheme. Rather than minimizing the squared op¬ 
erator distance D, we maximize the fidelity |4l| . which 
is in our context defined as an average overlap F = 

Tr(poutn). Here, a factorized system-bath input state 
Pin = !'!')(fkI (8 ) pb is transformed to gout = U by 
a total system-bath unitary operation U describing the 
device, after which an overlap is measured with a pro¬ 
jector n = f/p I'£')('!> I{7p. This projector is determined by 
the respective input system state |'I') transformed by the 
desired target system unitary Up and defines the mea¬ 
surement device. (• • ■) denotes the uniform average over 
all system states |'I'). Since the fidelity can thus be esti¬ 
mated from repeated binary measurements with sampled 
input states |'I') and gradually improved, this approach 
is independent of any assumptions about or knowledge 
of the environment and its coupling to the system. 

As an illustration, we consider again the task of im¬ 
plementing a Toffoli gate as shown in Fig. [6] In con¬ 
trast to Fig. [21 the three qubits j = 1,2,3 (again with 
randomly given Hj and Sj) are now coupled to a com¬ 
mon random environment (with randomly given Hb and 
Bj) rather than the described bosonic bath. Instead of 
Zermelo navigation, subsets of the qubits are brought in 
contact with the environment for fixed time periods = 
T in a fixed order given by the asynchronous protocol de¬ 
scribed above. The controls are the system-environment 
interaction strengths (cf. the above) over these fixed 
time periods. (Equivalently, we may fix the system- 
environment interaction strength and tune the time pe¬ 
riods Tk instead.) Plot 1 in Fig. [6] shows a simulation of 
the gradual increase of the fidelity as it would be seen in 


a control loop. For simplicity, we have assumed that dur- 



FIG. 6: Simulation of a closed loop gradient ascent of the 

-^;— —{V) 

fidelity F = Tr[l7(|v!')(v!'| 0 ^Bjtftll] , which can be esti¬ 
mated experimentally from sampled binary measurements of 
n = (vI'IHp, with the number of iterations i for three 

qubits sharing a random bath. Plot (1): Uf = Toffoli gate 
acting on the whole three-qubit space, 'P=T, Plot (2): Uf = 
Hadamard gate acting on a dynamically encoded decoherence 
free subspace V = span{|000), |111)}, using an asynchronous 
protocol with n = 576 time steps as controls. 

ing each iteration i, a reliable estimate of the fidelity has 
been obtained by a sufficient number of measurements. 
Depending on the respective environment and the num¬ 
ber of controls (time steps), such a control loop cannot 
guarantee to implement a desired gate in practice, as is 
illustrated in plot I of Fig. [HI for the Toffoli gate. Here, 
we have restricted ourselves to a 3-level ancillary system 
initially in a maximally mixed state pb = \l serving as 
environment, whose properties are unknown to the con¬ 
trol loop. Since with an idealized Markovian bath (i.e., a 
bath in a narrower sense) no change of the fidelity would 
be observable by definition, such a control loop can at 
the same time serve to measure - in terms of the change 
of F - the deviation of the environment from an idealized 
bath and thus to witness the presence of non-Markovian 
properties. The latter application does not require the 
realization of a unitary with perfect (unit) fidelity. 


B. Encoding a minimal noise subspace 

To mitigate or remedy this controllability problem, we 
may restrict ourselves to input states \^) belonging 

—i 

to a given subspace V. The fidelity F = Tr(poutn) is 

-(T) 

now given by the uniform average (• • •) with all states 
[fk) G F only. V thus adopts the role of a minimal noise 
or decoherence free subspace [d^, 1 ^ , which is dynami¬ 
cally encoded by the control loop, without knowledge of 










any symmetries in the open system dynamics. Plot 2 in 
Fig. ini illustrates this in the example of a Hadamard gate 
encoded in 'P=span{|000), |111)}. Due to the absence of 
any symmetries, the improvement of achievable fidelity is 
here due to the redundancy introduced by reducing the 
dimension of the (logical) quantum channel from 8 to 2. 

The asynchronous protocol is here used merely in con¬ 
nection to Sec. uni where the piecewise constant are a 
consequence of the concatenated Zermelo navigations. It 
is obvious that instead of this, three independent contin¬ 
uous time-dependencies Hj(t) can be applied just as well 
for the maximization of the fidelity, cf. e.g. |^ . 

V. SUMMARY AND OUTLOOK 

In summary, we have applied the solution for the quan¬ 
tum Zermelo navigation problem to a scenario where a 
set of open systems share a bath. We have shown that 
Zermelo navigation allows to generalize an analytic case 
of system-bath factorization of the total unitary evolu¬ 
tion operator to the presence of arbitrary background 
fields i?o, where only the individual system parts Sj in 
the coupling operators have to be navigated in their re¬ 
spective local winds Hj, whereas all bath operators re¬ 
main unaltered. We have given a quantum state inde¬ 
pendent refactorization time limit in terms of the mini¬ 
mal bath energy. Finally we addressed the reachability 
of general target unitaries by concatenations of Zermelo 
navigations, making use of the two-exponential form of 
the individual Zermelo unitaries, and demonstrated the 
feasibility in numerical examples. In addition, to ad¬ 
dress the navigation to a general target in an unknown 
environment, we have considered a measurement-based 
closed loop scheme. 

Interesting open questions we did not address and left 
for future work are the possibility of a symbolic solu¬ 
tion for quantum Zermelo navigation in a time-dependent 
wind Ho{t) that generalizes ([2]). From the numerical 
point of view, the gradient-based analyses in Secs. IIIII 
and HV] are thought as illustrations of proof of concept 
only (not addressing the question whether a local search 
of a target unitary starting at the identity yields a time- 


optimal navigation to this target unitary). For a com¬ 
prehensive investigation it would be interesting to incor¬ 
porate refined approaches such as subspace-selective self- 
adaptive differential evolution [s^ , Pareto front tracking 
[ 4 ^, two-stage hybrid optimization or an algebraic 
construction of the target unitary [ 43 . From the ex¬ 
perimental point of view, devising a realization of the 
modulation Q along with controlled interactions with 
an environment that obeys the desired bath properties 
is required. Potential candidates include Bose-Einstein 
condensates [l^, trapped ions [i^, atoms in an optical 
lattice d^, or arrays of optical cavities as depicted in 
Fig.m 

(a) (b) CA Csn-i 



FIG. 7: Possible schemes for controlled bath coupling: (a) 
static setup in which the systems (circles) are arranged in 
an array confined by a cavity U; (b) dynamic setup with 
two counterpropagating sequences of n — N/2 systems p (cir¬ 
cles) and q (dots) which interact during coincident passages 
through an array of 2n — 1 cavities: for identical interac¬ 
tions Kjk within the cavities, the term SjKjkSk then 

reduces to a collective coupling ~ PQ^ where P 

= Pj 8-nd Q — are variables of the respective 

system sequences. One may think of two trains of light pulses 
counterpropagating through a fiber, that may be regionally 
doped, or [(c), not shown] a ring resonator allowing repeated 
interactions. 
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